Search for the 0+ pentaquark in the reactions 7p K^K^n and 7p K^K^p 
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The exclusive reactions 7p — > R'^K^n and 7p K'^K'^p have been studied in the photon energy 
range 1.6-3.8 GeV, searching for evidence of the exotic baryon 0^(1540) in the decays nK^ 
and e+ pK°. Data were collected with the CLAS detector at the Thomas Jefferson National 
Accelerator Facility. The integrated luminosity was about 70 pb~^. The reactions have been isolated 
by detecting the K^ and proton directly, the neutral kaon via its decay to Ks — > n'^Tr~ and the 
neutron or neutral kaon via the missing mass technique. The mass and width of known hyperons such 
as E^, E~ and A(1116) were used as a check of the mass determination accuracy and experimental 
resolution. Approximately 100,000 A*(1520)'s and 150,000 </>'s were observed in the K^'K^n and 
K'^K'^P final state respectively. No evidence for the O"*" pentaquark was found in the nK^ or 
pKs invariant mass spectra. Upper limits were set on the production cross section of the reaction 
7P Q'^K'^ as functions of center-of-mass angle, nK'^ and pKs masses. Combining the results of 
the two reactions, the 95% C.L. upper limit on the total cross section for a resonance peaked at 1540 
MeV was found to be 0.7 nb. Within most of the available theoretical models, this corresponds to 
an upper limit on the width, Tq+, ranging between 0.01 and 7 MeV. 



I. INTRODUCTION 

The possible existence of baryon states beyond the 
usual qqq configuration is of fundamental importance 
for the understanding of hadronic structure. QCD does 
not prohibit the existence of exotic states with differ- 
ent configurations such as qqqqq. In fact, measurements 
of nucleon structure functions from high energy lepton- 
nucleon experiments have shown for example that "sea" 
quarks {qq pairs) contribute significantly to the total mo- 
mentum of the nucleon. Indeed usual baryons are ad- 
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mixtures of the standard qqq configuration and of qqqqq, 
qqqg, etc. 

In the past, experimental searches focused on pen- 
taquarks, i.e. baryons with a minimal ogggg structure. 
In 1997, Diakonov and collaborators made definite 
predictions about the masses and widths of a decuplet of 
pentaquark states (the so-called "antidecuplet" ) in the 
framework of a chiral soliton model. The most intriguing 
aspect of such a multiplct is the presence of three states 
with exotic quantum numbers or a combination of quan- 
tum numbers not allowed for ordinary baryons: the 0+ 
with S'=-|-l, the S and S"*" with S = —2. In particu- 
lar the positive strangeness 0+ is not compatible with a 
qqq state, requiring a minimal quark configuration of the 
type uudds. The widths of the exotic pentaquarks were 
predicted by this model to be very narrow (10-15 MeV) 
implying that if such states exist they should be directly 
visible in invariant mass spectra without the need for a 
more sophisticated partial wave analysis. 

The first evidence of a candidate was reported in 
October 2002 by the LEPS Collaboration, based on the 
re-analysis of existing data Several other experi- 

mental groups followed H HIl E S i, i, E EI 113 
reporting evidence of a peak in the mass range 1521- 
1555 MeV. The observation of an isospin partner of the 
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FIG. 1: The missing mass spectrum for the reaction 7p — > 
K^X. Peaks correspond to ground and excited states of weU- 
known hyperons. 



e+ (the e++) was recently reported by the STAR Col- 
laboration [13 while the observation of a second pen- 
taquark state, the S with dsdsu structure, was re- 
ported by the NA49 Collaboration Ti]. The first evi- 
dence for an anti-charmed pentaquark, Oc, was reported 
by the HI Collaboration [l^- While pentaquark signals 
observed in each experiment suffered from low statistics, 
the observations in many different reactions using differ- 
ent probes (photons, electrons, protons, neutrinos) and 
targets (protons, neutrons, nuclei) supported the pen- 
taquark's existence. On the contrary, subsequent re- 
analysis of data collected in a different set of experiments 

El 111 El El Sill 111 111 illp 111 11^ found 

no evidence of pentaquarks casting doubts about their 
existence. 

The experimental evidences, both positive and neg- 
ative, were obtained from data previously collected for 
other purposes in many reaction channels and in diverse 
kinematic conditions, thus may involve different produc- 
tion mechanisms. As a result, direct comparisons of 
the different experiments are very difficult, preventing 
a definitive conclusion about the pentaquark's existence. 
A second generation of dedicated experiments [H |^ , 
optimized for the pentaquark search, was undertaken at 
the Department of Energy's Thomas Jefferson National 
Accelerator Facility. These experiments covered the few 
GeV region in photon energy, where most of the positive 
evidence was reported, and collected at least an order of 
magnitude more data than used in the previous measure- 
ments. 

The exclusive reactions jp — > K^K^n and —> 
K"K°p were studied with the CLAS detector js^ with 



1.6 to 3.8 GeV energy photons, to look for evidence of 
the reaction jp K'^Q'^, where the decays into K'^p 
or K^n. The main results for the 7p K^K^n chan- 
nel were reported in Ref. [s^. In this paper we discuss 
in detail the analysis procedure and the results for both 
decay modes, combining them to give a final consistent 
result. 

For the 0+ nK^ decay mode, the measure- 
ment of all participating particles allows one to tag the 
strangeness of the reaction which clearly identifies the 
exotic nature of the baryon produced in association with 
the K'^ . For the other possible decay mode, 8+ pK'^, 
since a Ks was measured, the strangeness of the pKs 
invariant mass system is not defined. Nevertheless, we 
were motivated to analyze this channel since the major- 

^^^^f ^^^^Ti 

ity of the positive results 3, 8, 9, 10, 11, 12"| were reported 
looking at this decay mode. Moreover the CLAS accep- 
tance for the two reactions is somewhat complementary 
and the combination of the two channels results in com- 
plete kinematic coverage. According to many theoreti- 
cal predictions, e.g. IHI^IHIl^, the photon energy 
range covered by this experiment should be the best to 
explore since the 8"*" is expected to have its maximum 
cross section near the production threshold. Also, in this 
kinematic region, the CLAS detector provides a mass res- 
olution of few MeV and an accuracy in the mass deter- 
mination of 1-2 MeV, which is necessary to pin down the 
mass and width of any narrow peak in the spectrum. 

The jp K^K^n channel was previously investi- 
gated at ELSA by the SAPHIR collaboration in a 
similar photon energy range, finding positive evidence 
for a narrow 6+ state with M = 1540 MeV and FWHM 
r less than 25 MeV. The most recent analysis resulted 
in a total production cross section of the order of 50 nb. 
Since this experiment completely overlaps the kinematic 
regions of the SAPHIR experiment, the new results put 
those previous findings to a direct test for the first time. 
Results on pentaquark searches in the exclusive reaction 
7P — ^ K^K^p have never been published before. 

In the following, some details are given on the experi- 
ment (Secnj and its analysis fSec. lIIII) . The findings of 
7P K^K^n channel are compared with the SAPHIR 
experiment in Sec. II VI and the systematic checks are dis- 
cussed in Sec. [3 The final results are reported in the 
last section. 



II. THE EXPERIMENT 

The measurement was performed using the CLAS de- 
tector in Hall B at Jefferson Lab with a bremsstrahlung 
photon beam produced by a continuous 60 nA electron 
beam of Eq = 4.02 GeV impinging on a gold foil of thick- 
ness 8 X 10~^ radiation lengths. A bremsstrahlung tag- 
ging system |33 with a photon energy resolution of 0.1% 
Eq was used to tag photons in the energy range from 1.58 
GeV (about the 0^^(1540) production threshold) to a 
maximum energy of 3.8 GeV. A cylindrical liquid hydro- 
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FIG. 2: Hyperon production in the 7p R'^ K'^n reaction. 
Top: K'^ missing mass spectrum with the A* (1520). Bottom: 
nTr"*" (left) and nyr" (right) invariant mass spectra. The mass 
position and width of the measured peaks are shown (for the 
A* the experimental resolution is reported). For comparison 
the arrows show the world data value for the mass position. 
The shaded area indicates the retained events. 
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FIG. 3: Mass plots for analysis of the reaction 7p — > K^K'^n. 
Top: tt^tt" invariant meiss and the -K''^ peak. Bottom: missing 
mass for the reaction 7p K^K'^X and the neutron peak. 
The labels show the mass position and width of the measured 
peaks. For comparison the arrows show the world averages for 
the mass positions. The shaded area indicates the retained 
events. 



gen target cell 4 cm in diameter and 40 cm long was used. 
Outgoing hadrons were detected in the CLAS spectrom- 
eter. Momentum information for charged particles was 
obtained via tracking through three regions of multi-wire 
drift chambers 38] in conjunction with a toroidal mag- 
netic field (~ 0.5 T) generated by six superconducting 
coils. The polarity of the field was set to bend the positive 
particles away from the beam into the acceptance of the 
detector. Time-of-flight scintillators (TOF) were used for 
charged hadron identification j39|. The interaction time 
between the incoming photon and the target was mea- 
sured by the start counter (ST) This is made of 
24 strips of 2.2 mm thick plastic scintillator surround- 
ing the hydrogen cell with a single-ended PMT-based 
read-out. A time resolution of ~300 ps was achieved. 
The CLAS momentum resolution, ap/p is from 0.5 to 
1%, depending on the kinematics. CLAS is well suited 
for simultaneous multi-hadron detection as required by 
experiments searching for pentaquarks (this experiment 
required at least 3 hadrons detected). The detector geo- 
metrical acceptance for each positive particle in the rel- 
evant kinematic region is about 40%. It is somewhat 
less for low energy negative hadrons, which can be lost 
at forward angles because their paths are bent toward 
the beam line and out of the acceptance by the toroidal 
field. Coincidences between the photon tagger and the 
CLAS detector triggered the recording of the events. The 
trigger in CLAS required a coincidence between the TOF 
and the ST in at least two sectors, in order to select reac- 
tions with at least two charged particles in the final state. 
The collected data sample contains events from several 



reaction channels in addition to the reactions of interest. 
Reactions such as 7p —^ piT~^'!T~ , — > pto, jp — > K'^X, 
and jp have been used for system- 

atic studies. An integrated luminosity of 70 pb^^ was 
accumulated in 50 days of running in 2004. 



III. DATA ANALYSIS 

The raw data were passed through the standard CLAS 
reconstruction software to determine the 4-momenta of 
detected particles. In this phase of the analysis, cor- 
rections were applied to account for the energy loss of 
charged particles in the target and surrounding materi- 
als, misalignments of the drift chambers' position, and 
uncertainties in the value of the toroidal magnetic field. 
The energy calibration of the Hall-B tagger system was 
performed both by a direct measurement of the e~^e~ 
pairs produced by the incoming photons by ap- 

plying an over-constrained kinematic fit to the reaction 
7P — *■ p7r+7r~, where all particles in the final state were 
detected in CLAS 42]. The quality of the calibrations 
was checked by looking at the mass of known particles 
as well as their dependence on other kinematic variables 
(photon energy, detected particle momenta and angles). 
As an example, Fig. ^ shows the missing mass spec- 
trum of the reaction jp X: peaks of known hy- 
perons such as the A(1116), the S°(1193), and related 
excited states are clearly visible. 
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7P — > K^K^n after all cuts, obtained by integrating over all 
measured photon energies and K° angles. 



A. Pentaquark Analysis Strategy 

The data set was independently analyzed by three 
groups that made use of different analysis tools and pro- 
cedures. This strategy was adopted both to have a cor- 
roboration of the analysis results and an estimate of the 
systematic errors associated with the choice of the anal- 
ysis procedure. All three analyses agreed on the main 
conclusions. In the following sections we report in detail 
the analysis procedure of one group while the compari- 
son of the results from the different groups is discussed 
in Section rvOl 



B. Reaction identification: 7p K'^K^n 

The reaction 7p K^K'^n was isolated as follows. 
The K'^ was directly detected by the spectrometer, while 
the Ks component of the K'^ was reconstructed from the 
invariant mass of its tt+tt" decay (B.R. ~ 69%). The 
neutron was then reconstructed from the missing mass 
of all the detected particles. After its identification, the 
neutron mass used to calculate other kinematic variables 
was kept fixed at its nominal PDG value ^ ±3cr 

cut around the Kg and the neutron peaks was applied to 
isolate the exclusive reaction. A total of 320,000 events 
was selected by this procedure. Three background re- 
actions having the same final state as the reaction of 
interest were clearly identified: jp i4r+A*(1520) with 
A* (1520) ^ nK", 7p ^ X+S+tt", and jp -> K+j:~Tr+ 
with E^'^^^ ri7r+'^-'. Figure |2 shows the background 
hyperon peaks: the A* (1520) in the A'+ missing mass 
spectrum and the and peaks in the n7r+ and m:^ 



invariant mass spectra respectively. We found Mj2+ = 
1190±1 MeV and M^- = 1198±I MeV, with a measured 
experimental width a ~3.5 MeV for both of them. These 
are in excellent agreement with the world data and 
are a measure of the quality of the mass determination. 
Since these states have a much smaller width than the 
CLAS resolution, their observed widths provide an esti- 
mate of the experimental resolution. The reported values 
are in good agreement with the CLAS resolution esti- 
mated from simulations. To remove the contribution of 
these channels, a ±3(7 cut was applied around the S peaks 
while a ±24 MeV cut was appHed around the A* (1520) 
peak, resulting in a total of 160,000 retained events. The 
resulting Ks and neutron mass plots are shown in Fig.O 
where the two peaks are seen above small background 
levels. The peak positions of the reconstructed Ks and 
neutron masses were found to be 498±1 MeV and 939±1 
MeV respectively. 
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After all cuts, the resulting nK~^ invariant mass spec- 
trum is shown in Fig. ^ The spectrum is smooth and 
structureless. In particular no evidence for a peak or an 
enhancement is observed at masses close to 1540 MeV, 
where signals associated with the 0"*" were previously re- 
ported. To enhance a possible resonance signal not visi- 
ble in the integrated spectra, we assumed the two-body 
reaction 7p A'"8"''(1540) and selected different K^ 
center-of-mass angle intervals. Fig. |31 shows the nK^ 
invariant mass spectrum for different cos^^J^ ranges. 
Monte Carlo studies of the CLAS acceptance for this re- 
action (see Sec. nil E|) showed that we could detect events 
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over the entire angular range, in spite of a drop in the ef- 
ficiency at forward K'^ angles, O'^g^ < 30°. No structures 
were found in the spectra when specific angular ranges 
were selected. 

As a demonstration of our sensitivity to baryon reso- 
nances, we derived the A* yield. We fit the miss- 
ing mass spectrum of Fig. |5] by a Breit-Wigner function 
convoluted with a Gaussian function to account for the 
detector resolution, plus a second-order polynomial back- 
ground. To derive the A* yield, following Ref. (4^, we 
used the Breit-Wigner form: 



BW = 



Momr(m) 



(m2 - M2)2 + M2r(m)2 ' 



(1) 



where: 



r(m) = To I 



21+1 



Mq and Fq are the resonance mass and intrinsic width, Q 
is the momentum in the rest frame of the n — K'^ sys- 
tem, Qo is the same quantity evaluated at the A* (1520) 
peak, and / is the n — relative orbital angular mo- 
mentum [1 — 2 for the A*). In the fit Fq was fixed to 
15.6 MeV ^3 while the a of the Gaussian function was 
allowed to vary. Integrating the Breit-Wigner line in the 
mass range 1.45-2.0 GeV we obtained a A*(1520) yield of 
99, 000 ± 10, 000. The quoted error is dominated by the 
systematic uncertainty related to the shape of the Breit- 
Wigner and the underlying background. The mass posi- 
tion was found to be 1518 ± 2 MeV, in good agreement 



with world data |43l |. while the experimental resolution 
was found to be ~ 5 MeV, typical for CLAS ^] . 



Reaction identification: 



IP 



K''K°p 



To isolate the reaction jp K'^K^p, the proton was 
detected by the spectrometer and a Ks meson was re- 
constructed from the invariant mass of its tt+tt^ decay 
(B.R. ~ 69%). The second neutral kaon was then re- 
constructed from the missing mass of all the detected 
particles. A ±3(T cut around the Ks and the missing 
kaon peaks was applied to isolate the exclusive reaction. 
A total of 750,000 events were selected by this proce- 
dure. The reaction ^p — > p(j) with (p —f K^Ks has the 
same final state as the reaction of interest. The reaction 
7P A(1116)7r+iir° "> pTT^Tr+if^ also contributes to the 
background. Fig. El shows the background peaks: the 4> 
shows up in the proton missing mass spectrum and the 
A(1116) in the p7r~ invariant mass spectrum. We found 
= 1019 ± 1 MeV and Ma = 1116 ± 1 MeV, with a 
measured experimental resolution a of about 7 MeV and 
2 MeV respectively. The cj) peak was fitted with a Breit- 
Wigner (Fq — 4.2 MeV), convoluted with a Gaussian 
describing the CLAS resolution, in the same way as was 
done for the A* (1520). The obtained masses are in agree- 
ment with world data [4^ . To remove the contribution 
of these channels, only events with a Mxijp ^ pX) > 
1.04 GeV and M(pTT~) > 1.13 GeV were retained, result- 
ing in a total of 550,000 events. The two neutral kaon 
mass spectra after the background rejection are shown 
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FIG. 8: The pKs invariant mass (top) and Ks missing mass 
(bottom) for the reaction 7p —> K^K'^p after all cuts, ob- 
tained by integrating over all measured photon energy and 
angles. Shoulders corresponding to excited S states are visi- 
ble in both spectra. 



in Fig. The two masses were found to be respectively 
498±1 MeV and 496±3 MeV. 

The selected final state contains two neutral kaons, one 
detected and one missing, therefore a possible G"*" peak 
can show up in two ways: in the invariant mass spectrum 
of the pKs system or in the missing mass spectrum of the 
detected Ks- Fig. [Slshows the two spectra after all cuts: 
both of them are smooth and structureless. In particular 
no evidence for a peak or an enhancement is observed at 
masses close to 1540 MeV, where signals associated with 
the @^ were previously reported. To enhance a possible 
resonance signal not visible in the integrated spectra, we 
assumed the two body reaction 773 A'°8+(1540) and 
selected different K'^ center-of-mass angle intervals. The 
angle was calculated using the reconstructed kine- 
matic variables of the missing kaon in the first case and 
the measured kinematic variables of the observed kaon 
in the second case. Fig. shows the pKs invariant mass 
and the Kg missing mass spectra for forward and back- 
ward cos 9'^^ ranges separately. When the pKs system is 
considered, Monte Carlo studies showed that the CLAS 
acceptance is maximum at forward 9^^^ and therefore 
complementary to what we found for the — *■ nK'^ 

No structures were 



decay mode (as seen in Sec. IIII B|) 
found in any of the spectra when specific angular ranges 
were selected. 

As shown in Fig. [T] the A'5 peak sits over a large 
background mainly related to multi-pion production. A 
cleaner sample is obtained by applying a cut on the Kg 
decay length: in fact, due to the sizeable Ks mean life 
(cT ~ 2.68 cm), its decay vertex {Ks ^ tt+tt^) is dc- 
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FIG. 9: The pKs invariant mass (left) and Ks missing mass 
spectra (right) for the reaction 7p K'K'^p for different an- 
gular ranges: cos ^^2^ < (top), and cos 0^5^ > (bottom). 
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FIG. 10; Effect of the Ks decay length cut on the analysis. 
Top: tt^tt" invariant mass and the Ks peak (as in Fig|7Jtop). 
Bottom: pKs invariant mass (as in FiglHJtop). 



tached from the primary production vertex {pK'^K'^). 
Taking into account the vertex resolution of the CLAS 
detector (~ 0.3 cm) and the Ks ct, we applied a 3 cm 
decay length cut, obtaining the mass spectra of Fig. 1101 
Despite the use of a cleaner Ks sample (top panel), no 
structures are present in the pKs invariant mass spec- 
trum (bottom panel), confirming the results reported 
above. The Ks decay length cut improves the signal- 
to-background ratio for the Ks identification, cleaning 
the data sample from the multipion contamination. On 
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11: The nK'^ invariant mass spectrum for the reaction 
7P — > K°K^n fitted for a 9+ mass of 1540 MeV using the 
three different methods described in the text: from top to 
bottom the fit of the histograms, the likelihood fit of the 
unbinned events, and the background fit with the exclusion 
of the Q'^ mass window (indicated as hatched area). The 
obtained yields for signal and background are shown on top 
of each plot. Insets show the zoom of the mass region 1520- 
1560 MeV. 



the other side it reduces the Ks yields by almost a factor 
five and, due to the strong correlation of the Ks mean 
life with the kaon momentum, it results in a momentum- 
dependent cut on the Ks sample, difficult to reproduce 
by Monte Carlo simulation. This also distorts the pKs 
invariant mass spectrum at low values where a possible 
resonance is more likely produced. For these reasons, the 
upper limit for this decay mode are evaluated without the 
vertex cut. 



D. Upper limit on the O"'" yields 

Since no signal was observed, an upper limit on the 0+ 
yield was extracted and transformed to an upper limit on 
the production cross section in each reaction channel. In 
this section we discuss in detail the procedure adopted 
for the channel —> K^K'^n and, since the procedure 
is the same, we summarize the results for the channel 
7P K°K°p. 

The was assumed to be a narrow peak over a 
smooth background. The nK~^ invariant mass spectrum 
was fit to the sum of a Gaussian-shape resonance and a 
fifth-order polynomial representing the background. The 
resonance position was varied from 1520 to 1600 MeV in 
5 MeV steps while the width was kept fixed assuming 
the dominance of the experimental resolution over the 
intrinsic width as suggested from recent analyses of KN 
scattering data 45]. 

The 0+ width, aQ+ , was derived by means of a Monte 
Carlo simulation. A zero-width resonance was gener- 
ated and projected over the CLAS detector, applying 
the same analysis chain used to process the data. A 
width of (Tq+ of 3-4 MeV was obtained, which weakly de- 
pends on the photon energy and the i^" emission angle. 
To check the consistency of the experimental resolution 
obtained from the Monte Carlo simulations, the same 
procedure was applied to the reactions 7p if+E+7r~, 
and 7p — > A'^^S^tt"*" where a direct comparison to the 

The width of the two 



data is possible (see Sec. IIII B|) 
hyperons derived from the simulations was found to be 
compatible with the measured values. For each value of 
the 0^ mass, the resonance and background yields were 
extracted using three diS'erent fit procedures: (1) a 
fit of the mass distribution binned in 4 MeV channels; 

(2) a likelihood fit of the unbinned nK^ spectrum; (3) 
as in case (1) but with the background function being 
fit after excluding the signal region defined as ±3(Te+ 
around the selected mass value. In all cases the back- 
ground yield was obtained by integrating the polynomial 
function over ±3(Te+ of the selected 0+ mass value. In 
methods (1) and (2) above, the signal yield was obtained 
as the integral of the resulting Gaussian, while in method 

(3) it was obtained as the difference between the num- 
ber of observed events and the background integrated 
over ±3cre+ of the chosen 0+ mass value. The same 
procedure was then repeated subdividing the data into 
16 cos 0^2^ bins producing binned spectra. The results 
of the three methods applied to the integrated spectrum 
are shown in Fig. ^2 In general, the signal yields ob- 
tained with the three procedures are compatible with 
zero within 1 or 2 sigma, confirming that no evidence 
for 0+ production is observed in the mass range 1.52-1.6 
GeV. The results of the binned fi^ts and un-binned like- 
lihood fit are in good agreement with each other showing 
that the binning effects are small. The measured yields 
and the background are shown in the top and middle pan- 
els of Fig. El In general the results of the three methods 
(x^, likelihood, and without the 0+ mass window) 
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FIG. 12: Signal (top) and background (middle) yields ob- 
tained with the three fit procedures for the reaction 7p 
K^K^n: fit of the histograms (circles), likelihood fit of 
the unbinned events (squares), and background fit with the 
exclusion of the mass window (triangles). The dashed 
lines are to guide the eye. The three results were combined 
together taking the average of the event yields (solid line) 
and transformed in 95% C.L. upper limit shown in the bot- 
tom panel. 



are consistent with each other as expected by the domi- 
nance of the background over the signal yield. They were 
combined by taking the average of the event yields, for 
both signal and background, in the conservative assump- 
tion of totally correlated measurements. The averaged 
yields were transformed into upper limits of the true 
yield using the Feldman and Cousins approach . This 
method determines proper confidence level boundaries 
for small signals over a background taking into account 
external constraints (e.g. the true yield is constrained to 
be positive). In addition it decouples the goodness-of-fit 
confidence level from the confidence level interval. The 
resulting upper limit at 95% C.L. is shown in the bottom 
panel of Fig. 1121 it is almost flat with the maximum value 
around 1545 MeV. 

On average, the upper limit at 95% C.L. on the yields is 
A^e+ ~ 220. The ratio of the yield of 9+ to the A* (1520) 
has also been obtained: j;t^ 220/99000 - 0.22% 

™A*(1520) 

(the two yields are not corrected for the CLAS efficiency). 
The same procedure was repeated for the 7p K^K^p 
channel to derive an upper limit at 95% C.L. on the yield 
as a function of the pKs invariant mass and on the dif- 



FIG. 13: The same as in Fig.[Tlfor the 7p K°K°p channel. 
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FIG. 14: Binned CLAS efficiency for the reaction 7p 
^""0^. Top: nK^ assuming the exchange domi- 

nance hypothesis; bottom: O"'" pK'^ assuming the calcula- 
tion of Ref. 34] with no is'* -exchange. These two correspond 
to the 'worse case scenario' for each decay mode respectively. 

ferential cross section assuming a 0+ mass of 1540 MeV. 
The upper limit was derived from the pKs mass spectrum 
only (the Ks missing mass spectrum was ignored because 
it is correlated to the pKs mass spectrum being built 
from the same event sample). The better CLAS accep- 
tance for the proton and the Ks coming from the 0+ de- 
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FIG. 15: The 95% Confidence-Level upper limit on the total 
cross section for the reaction jp K'^Q'^ with 0^ nK'^ 
(top) and 0"*" pK^ (bottom). The dashed line is to guide 
the eye. 



cay leads to a complete angular coverage complementary 
to the 7p — > K'^K'^n channel. The O"*" was searched for 
as a narrow peak by fitting the pKs mass spectrum with 
a Gaussian curve with (Tq+ ~4 MeV inferred from ded- 
icated Monte Carlo simulations, plus a fifth-order poly- 
nomial representing a smooth background. The mass re- 
gion 1520-1600 MeV was scanned in 5 MeV steps. To 
derive the yields of a possible resonance and associated 
background the three fit procedures described for the 
7P K^K'^n channel were applied to the integrated 
spectrum and the cos 0^J^-binned spectra for a fixed 6 
mass of 1540 MeV. The three results were combined tak- 
ing the average of the event yields. These were trans- 
formed into 95% C.L. upper limits. Figure [T^ shows the 
comparison of the 3 fitted yields as a function of the pKs 
invariant mass and the resulting 95% C.L. upper limit. 



E. Upper limits on the 0+ Cross Section 

The 95% C.L. upper limits on the yield described in the 
previous section were then transformed into limits on the 
production cross section according to the following 
formula: 



K+ 



C 



'lK+ 



pKO 



(2) 



(3) 



where N is the 95% C.L. limit on the 8+ yield, e is 
the CLAS efhciency, £ is the integrated luminosity, h 
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FIG. 16: The 95% Confidence-Level upper limits on 
the differential cross section da / d cos 9'^o' for the reaction 
jp _R'''0+(154O) obtained from the two decay modes: 
0+(154O) nK+ (top) and 0+(154O) ^ pK" (bottom). 
The dashed line is to guide the eye. 



is the branching ratio for the O"*" decay, and subscripts 
indicate the decay mode of the Q^. Each branching ratio 
is assumed to be 50%. The branching ratios for neutral 
kaon decay to K° ^ Ks ^ tt+tt" (50% • 68%) were 
included in the CLAS efficiency. 

The luminosity was obtained as the product of the 
target density and length, and the incoming photon 
flux measured during the experiment corrected for data- 
acquisition dead time. When the mass was varied 
from 1.5 to 1.6 GeV, the production threshold in beam 
energy moved from 1.65 to 1.85 GeV. The photon flux 
used in the cross section estimate was calculated accord- 
ingly. 

The CLAS detection efficiency was obtained by means 
of detailed Monte Carlo studies. The reaction 
^°9+ and subsequent G"*" decay to nK~^ and pK^ was 
generated assuming different production mechanisms: 
t— exchange dominance (the K'^ mainly produced at for- 
ward angles in the center-of-mass system), u— exchange 
dominance {K^ at backward angles), cos^^^ uniformly 
distributed, and using the predictions of the model in 
Ref. [s^. For the i— exchange hypothesis we assumed 
the same angular distribution as for the reaction 7p 
A*(1520)fir"'", which exhibits the typical i— channel for- 
ward peaking behavior (approximately an exponential 
with a slope of —2.5 GeV~^). The u— exchange distribu- 
tion was generated the same way except that the center- 
of-mass angles of the K'^ and 0+ were interchanged. 

For 7p K°e+, 6+ K+n, the CLAS overall de- 
tection efficiencies obtained with the various production 
mechanisms vary between 1% for the exchange hy- 
pothesis to 1.8% for the angular distribution of Ref. |3J| 
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FIG. 17: The 95% Confidence-Level upper limits on the to- 
tal cross section as a function of E.y for the reaction 7p — > 
i^°e+(1540) with e+ nK+. The dashed line is to guide 
the eye. 
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FIG. 18: The CLAS missing mass showing the 

A* (1520) peak before (solid) and after (dashed) cutting on 
cos6i|f >0.5. 



when no K* exchange process is included. As a func- 
tion of 9'^^^ all the different hypotheses gave a com- 
parable efficiency: almost flat from 180° to 90° (about 
2%) and then smoothly dropping at forward angles. For 
7P X^G+j 8+ — > the efficiency varied between 

1% for the angular distribution of Ref. [s^] with no K* 
exchange process to 1.8% for the i— channel hypothesis, 
with an angular dependence complementary to the other 
channel (smoothly increasing from backward to forward 
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FIG. 19: The CLAS nK'^ invariant mass with all the cuts 
as in Ref. 5]. The inset show the same spectrum with the 
binning used by the authors of the same reference. 



angles) . For each branch, the model that yielded the low- 
est efficiency was chosen for conservatism. The resulting 
efficiencies are shown in Fig. ^1 

The upper limits on the total cross sections as a func- 
tion of the 0^ mass were obtained independently for 
each decay mode using the model assumptions described 
above. The results are shown in Fig.^l For 0+ nK^ 
a 95% C.L. upper limit of 1.0 nb was found for Mq+ = 
1540 MeV. The corresponding limit for 0+ pK^ was 
1.3 nb. 

The 95% C.L. upper hmit on the 6+ (1540) differen- 
tial cross section da /d cos is shown Fig. ^1 using 
the same assumption on the production mechanisms as 
for the evaluation of the upper limit on the total cross 
section. However, for this quantity no significant differ- 
ence was found when the other hypotheses were used in 
the efficiency evaluation. For the nK'^ decay mode the 
cross section limit remains within 1-2 nb for most of the 
angular range and rises at forward angle due to the re- 
duced CLAS acceptance. For the pK° decay mode, the 
cross section limit is within 2-5 nb over the entire angular 
range. 

Finally, for the reaction jp /i°6+(1540) with 6+ 
nK^, the 0"'"(154O) total cross section upper limit as a 
function of the photon energy is shown in Fig. 1171 the 
behavior reflects the CLAS efficiency, which is reduced 
at low energy near the 8"*" (1540) production threshold, 
and then increases with energy, resulting in a better limit 
for higher energies. 
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IV. COMPARISON WITH EXISTING DATA 

Our result for the reaction 7p — > K^Q~^ K'^K^n is 
in clear disagreement with the findings of Ref. which 
reported a signal of 55 events at a mass of 1540 MeV 
corresponding to an estimated total cross section of 50 
nb. 

In order to better compare our data with this exper- 
iment, the kinematic cuts described in Ref. @ were ap- 
plied to the present data. The photon energy was limited 
to be below 2.6 GeV and a possible forward peaked pro- 
duction of the 0+ was selected by applying the angular 
cut cos 6*^0^ > 0.5. This cut also reduces the hyperon 
production yield. The effect is shown in Fig. 1181 where 
the missing mass of the is shown before and after 
the angular cut: the A* (1520) is clearly suppressed. As 
a result, no hyperon rejection cuts were applied. 

The K'^ missing mass spectrum, after all cuts, is shown 
in Fig. ^1 with two different bin sizes reflecting the CLAS 
and the SAPHIR resolutions. There is no evidence of a 
6+ peak. 

Applying the same procedure described above, we eval- 
uated a 95% confidence level limit on the 8+ (1540) yield 
with SAPHIR selection cuts of 90 events. 

To derive the A(1520) yield, the K'^ missing mass spec- 
trum obtained before the angular cut was fit by a Breit- 
Wigner function plus a second-order polynomial back- 
ground, with the same procedure described in Sec. 1111 lj| 
obtaining (57,000±5,500) A*(1520)s. This number is to 
be compared to 630±90 reported by SAPHIR Ref. Q. 
The ratio between observed 0+ and A(1520) in this ex- 
periment is ^ 0.16% differing by more than a factor 50 
from the value quoted in Ref. 5] . All the yields reported 
above are notcorrected for detector acceptances. 

No results were published on the search of the 0+ in 
the exclusive reaction 7p K'^K'^p and therefore no 
comparison is possible for this channel. 



V. SYSTEMATIC ERRORS 

In the evaluation of the upper limit we have considered 
the following sources of systematic errors: determination 
of the mass resolution for the 0"*" resonance, determina- 
tion of the signal and background yields from the mass 
spectra, evaluation of the CLAS efficiency, detector in- 
efficiencies, photon flux normalization, and dependence 
on the analysis procedure. The first three sources were 
already included in the quoted upper limit as explained 
in the following subsections. In particular, the model 
dependence in evaluating the CLAS efficiency was esti- 
mated by comparing the results obtained using different 
models for the production cross section. The resulting ef- 
ficiencies differ by a maximum of a factor two. The upper 
limits were estimated using always the worse case sce- 
nario. The remaining sources of systematic uncertainty, 
summarized in Table ^ result in an overall systematic 



error of e = 25% accounted for by multiplying the upper 
limit by (1+e). 

TABLE L Systematic errors on the upper limit evaluation. 



Source 


Error (%) 


Detector inefficiencies 
Photon flux normalization 
Analysis procedure 


10 
10 
<20 



In the following, the different contributions to the total 
systematic error are discussed in more detail. 



Mass resolution and evaluation of signal and 
background yields 



As discussed in Sec. IIII Dl the mass resolution for a 
narrow resonance was estimated from Monte Carlo sim- 
ulation; the reliability of simulations in reproducing the 
CLAS resolution was tested comparing the observed reso- 
lution for known narrow resonances and a maximum dis- 
crepancy of 20% was found. The resolution for the 0+ 
peak extracted from Monte Carlo was therefore rescaled 
by a factor 1.2 to account for this. 

The comparison of different fitting procedures provides 
an estimate of the associated systematics. The upper 
limits were derived combining the results of the three ffis, 
including therefore an estimate of the associated error. 
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FIG. 20: The final nK'^ spectra for the reaction jp 
K^n obtained by the three independent analyses. 



B. Detector inefficiencies and normalization 

As a check of the accuracy of the CLAS detector sim- 
ulations and photon flux normalization, the differential 
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FIG. 21: Signal (top) and background (bottom) yields ob- 
tained by the three independent analyses for the reaction 
7P K^K'^n. Different symbols refer to the different analy- 
ses. The solid line shows the average of the three results. 
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FIG. 22: The 95% C.L. upper limit on the cross section for the 
reactions 7p ^ K°Q+ K^K+n (top) and 7p ^ K°Q+ 
K^K'^P (bottom) derived by the combination of the three 
analyses. The dashed line is to guide the eye. As explained 
in the text, the three analyses were largely independent. 



and the total cross section for several known reactions 
were derived from this data set. Due to the high multi- 
plicity of charged particles in the final state in our data 
(similar to K'^K^n and K'-'K'^p) and the existence of 
precise measurements that can be taken as a reference, 
the reactions 7p — > puj and 7p — > K'^ A{1116) were used 
as a test of the different ingredients used in our analy- 
sis. Moreover, the measurement of different final states 
of the same reaction, such as p7r"'"(7r~7r*') and p7r"'"7r~(7r°) 
for the uj photoproduction and K~^p{'K~) and K^pi:^ 
for the iir+A(1116) channel, were also used to test the 
hardware trigger, the photon flux normalization and the 
procedure to extract the CLAS efficiency. The differen- 
tial and the total cross sections extracted from this data 
set agree with each other and with the world data within 
the experimental error, verifying that the different steps 
in the analysis are in control at the 15% level. 

To directly check our ability to observe the final state 
involved in the pentaquark decaying into the nK^ fi- 
nal state, the cross section for the reaction 7p 
i4'+A*(1520) K+K°n was also extracted. As shown 
above, the A* (1520) peak is clearly visible, and the sam- 
ple of 100,000 A*'s made possible an analysis deriving 
both the differential and the total cross sections. The 
A* (1520) cross sections we obtained were compared with 
data from the SAPHIR [13 and NINA ^ collabora- 
tions. We are in good agreement with the SAPHIR re- 
sults but find a much lower cross section than that re- 
ported by NINA. The results of all these measurements 
will be reported elsewhere. 



C. Dependence on the analysis procedure 

Three independent analyses were conducted both on 
the reactions of interest for the pentaquark search and 
on reference reactions as -fp A* {1520)K~^ K^K'^n. 
This enabled an evaluation of the systematic errors asso- 
ciated with the analysis procedures and provided a cross 
check on the results. The three analyses were not totally 
independent since they all used the same raw data and 
the same basic corrections to the measured kinematic 
quantities such as the energy loss. However, they used 
different particle identification schemes, different detec- 
tor calibration procedures (for both the tagger system 
and the CLAS spectrometer) , different Monte Carlo sim- 
ulations to evaluate the CLAS efficiency, and different fit 
procedures to extract yields. 

The A* (1520) differential and total cross sections ob- 
tained by the three analyses were found to be consis- 
tent with each other at the 10% level and agree with the 
SAPHIR measurement at the same level. 

For the reaction jp K'^Q'^, all three analyses agreed 
on the main conclusion: the nK^ and pKs spectra are 
smooth and structureless and, in particular, no signal 
is observed at 1540 MeV where the 8+ pentaquark has 
been widely reported. 

The nK'^ mass spectra for the reaction jp — + K^K^n 
obtained by the three analyses are shown in Fig. EOl The 
difference in the shape of the spectra is mainly related to 
the different particle identification schemes adopted by 
the three groups while the small differences in the bin- 
to-bin fluctuations are due to the different kinematic cor- 
rections applied in the analyses. Signal and background 
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yields as a function of the 0+ mass resulting from the fit 
of the three spectra are shown in Fig. |^ The three anal- 
yses are consistent at the 10% level on the background 
estimate while they differ in the event yield evaluation. 
This discrepancies, reflecting the different choices in the 
analysis procedures, provide an estimate of the system- 
atic error associated to the extraction of the upper limits. 

The same comparison was performed for the 7p — > 
K^K^p channel with similar results. 



VI. RESULTS 

A. Upper limits on the production cross section 

The independent analyses were combined, taking the 
average of the event yields, and transformed into the 95% 
C.L. upper limit on the yields with the Feldman and 
Cousins procedure. They were then transformed into the 
95% C.L. on the cross section using the CLAS efficiency 
evaluated in the most conservative scenario. Results are 
shown in Fig. 

For the reaction 773 — ^ A'"9+ with 8+ nK^ the 
upper limit at 95% C.L. on the 0+ production cross sec- 
tion varies between 0.5 nb and 1.3 nb as a function of 
the nK~^ invariant mass with a value of ^--^0.8 nb at 1540 
MeV. For the reaction jp K"Q+ with 9+ pK° 
the upper limit at 95% C.L. on the 0+ production cross 
section varies between 0.5 nb and 2.5 nb as a function of 
the pKs invariant mass with a value of ^1.5 nb at 1540 
MeV. The results for the two decay modes are similar in 
value and set stringent upper limits on the models which 
predict these long-lived pentaquark states. 

So far the two reaction channels were independently 
analyzed. Assuming they result from the two possible 
decay modes of the 0+, they can be combined to give a 
single upper limit. 

As shown in the previous sections, we estimated a sig- 
nal yield Si (with i = 1 and 2, corresponding to nK^ 
and pK'^ branching mode), e.g. the area of a Gaussian of 
fixed width and fixed mass fit to the mass histogram, and 
a background yield Bi, e.g. the polynomial background. 

These were corrected for the detection efficiency and 
luminosity to obtain the two corresponding cross sections 
((Ti) and associated errors {6i): 

Si \J Si + Bi 

Oi — , Oi — 

C ti hi C ei bi 

with hi the corresponding branching ratios, C the inte- 
grated luminosity and the CLAS detection efficiency. 

The cross section for a possible 0+ is then built as the 
weighted average of the two, using the CLAS efficiency 
evaluated using the five models described in Sec lIII El 
The dependence on the model assumptions resulted to 
be within 30%. The largest cross section was obtained 
in the hypothesis of exchange dominance then chosen 
for conservatism. 



TABLE II: The 6+ total cross section (nb) predicted by dif- 
ferent models assuming = 1 MeV. 





1/2+ 


1/2- 


3/2+ 


3/2" 


[34] 


100 


0.4 






f35] 


0.22 


0.1 






[36] 


2.0 




1.0 


3.0 


[49] 


6.9 


3.4 


3.2 


17.7 


[50] 


15.0 









TABLE III: Upper limits on the 9+ width (MeV) assuming 
a 95% C.L. of 0.7 nb on the 0+ total cross section. 





1/2+ 


1/2- 


3/2+ 


3/2" 


[34] 


0.01 


1.8 






[35] 


3.2 


7.0 






[36] 


0.35 




0.7 


0.23 


[49] 


0.10 


0.2 


0.2 


0.04 


[50j 


0.05 









B. Upper limits on 

The 0+ production cross section is directly connected 
to the 0+ width (see for example Ref. H^)- There- 
fore upper limits on the cross section imply upper lim- 
its on the resonance width. However, this connection 
depends strongly on the theoretical model, differing by 
more than an order of magnitude for the available calcu- 
lations 34, 35, 36, 49, 50] . In TablellTl we summarize var- 
ious theoretical predictions for the total cross section av- 
eraged in the energy range 2-4 GeV for different assump- 
tions for parity and spin of the 0+ and ~ 1 MeV. 
For example, assuming = 1/2+ our upper limit of 
0.7 nb on cr{-fp K°Q+) results in a re+ < 3.2 MeV 
within the Regge approach of Ref. and < 0.35 
MeV for the other models. The upper bounds on the 0+ 
width for these models are summarized in Table Hill 

VII. CONCLUSIONS 

In this paper we report the results of the first Jefferson 
Lab high statistics and high resolution experiments en- 
tirely devoted to the pentaquark search on a nucleon tar- 
get. The reactions 7p K'^K'^n and ^p K^K^p were 
studied in search of evidence of the 0+ pentaquark in 
the nK'^ and pK^ decay channels. The final states were 
isolated by detecting the or proton, the Ks via its 
decay to tt+tt^ and identifying the neutron or the second 
neutral kaon with the missing mass technique. For the 
former decay mode, the direct measurement of the al- 
lows the definition of the strangeness of any resonance ob- 
served in this final state. A total of 160,000 and 550,000 
events were selected for the reaction 7p K'^K^n and 
7P K^K^p respectively, after the exclusion of back- 
ground reactions. The 0+ was searched for as a nar- 
row resonance in the nA'+ and pKs mass spectra with 
a width of 3-4 MeV corresponding to the CLAS reso- 
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FIG. 23: The 95% C.L. upper limit for the reaction 7p 
K^O'^ combining the two channels jp K'^K'^n and jp 
K'^K^p. The 95% C.L. upper limit on the O"*" cross section 
as a function of O"*" mass (top) and the 95% C.L. upper limit 
on the differential cross section da / d cos O'^o' for a fixed 0^ 
mass of 1540 MeV. The dashed line is to guide the eye. 



lution for these channels in the kinematic region of this 
experiment. Both mass spectra were found to be smooth 
and structureless. No evidence for a narrow resonance 
was found in the mass range 1520-1600 MeV. Combin- 
ing the results of the two decay modes, we set an upper 
limit of 0.7 nb (95% confidence level) on the total pro- 
duction cross section for the reaction 7p A'°6+(1540). 
This contradicts the results previously reported for a res- 
onance in the reaction channel 7p — > K^K^n. 

The accuracy in the mass determination was found to 
be 1-2 MeV from the comparison of the measured masses 
of known particles with world data. The quality of the 
data and the analysis procedures were tested by deriving 
the differential and the total cross section for some known 
reactions and obtaining an agreement within the exper- 
imental errors with existing measurements. The same 
conclusions were found by several independent analyses, 
giving confidence in our final results. 

Due to the loose hardware trigger of the experiment, 
the same data set was analyzed to study the reactions 
7P ^ K*Q+ and 7p -> K-Q++ [5l|. These findings, 
together with the results coming from other pentaquark 
search experiments at Jefferson Lab, could clear up the 
debate about the existence of the pentaquark. 
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